Search for: All records

We present an adjoint-based optimization method to invert for stress and frictional parameters used in earthquake modeling. The forward problem is linear elastodynamics with nonlinear rate-and-state frictional faults. The misfit functional quantifies the difference between simulated and measured particle displacements or velocities at receiver locations. The misfit may include windowing or filtering operators. We derive the corresponding adjoint problem, which is linear elasticity with linearized rate-and-state friction and, for forward problems involving fault normal stress changes, nonzero fault opening, with time-dependent coefficients derived from the forward solution. The gradient of the misfit is efficiently computed by convolving forward and adjoint variables on the fault. The method thus extends the framework of full-waveform inversion to include frictional faults with rate-and-state friction. In addition, we present a space-time dual-consistent discretization of a dynamic rupture problem with a rough fault in antiplane shear, using high-order accurate summation-by-parts finite differences in combination with explicit Runge–Kutta time integration. The dual consistency of the discretization ensures that the discrete adjoint-based gradient is the exact gradient of the discrete misfit functional as well as a consistent approximation of the continuous gradient. Our theoretical results are corroborated by inversions with synthetic data. We anticipate that adjoint-based inversion of seismic and/or geodetic data will be a powerful tool for studying earthquake source processes; it can also be used to interpret laboratory friction experiments. 

Abstract Explosive volcanic eruptions radiate seismic waves as a consequence of pressure and shear traction changes within the conduit/chamber system. Kinematic source inversions utilize these waves to determine equivalent seismic force and moment tensor sources, but relation to eruptive processes is often ambiguous and nonunique. In this work, we provide an alternative, forward modeling approach to calculate moment tensor and force equivalents of a model of eruptive conduit flow and chamber depressurization. We explain the equivalence of two seismic force descriptions, the first in terms of traction changes on conduit/chamber walls, and the second in terms of changes in magma momentum, weight, and momentum transfer to the atmosphere. Eruption onset is marked by a downward seismic force, associated with loss of restraining shear tractions from fragmentation. This is followed by a much larger upward seismic force from upward drag of ascending magma and reduction of magma weight remaining in the conduit/chamber system. The static force is upward, arising from weight reduction. We calculate synthetic seismograms to examine the expression of eruptive processes at different receiver distances. Filtering these synthetics to the frequency band typically resolved by broadband seismometers produces waveforms similar to very long period seismic events observed in strombolian and vulcanian eruptions. However, filtering heavily distorts waveforms, accentuating processes in early, unsteady parts of eruptions and eliminating information about longer (ultra long period time scale depressurization and weight changes that dominate unfiltered seismograms. Our workflow can be utilized to directly and quantitatively connect eruption models with seismic observations. 

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.